1386 


IEEE  TRANSACTIONS  ON  COMPUTER-AIDED  DESIGN  OE  INTEGRATED  CIRCUITS  AND  SYSTEMS,  VOL.  19,  NO.  1 1,  NOVEMBER  2000 


Comments  on  “Sympathy:  Fast  Exact  Minimization  of  Fixed 
Polarity  Reed-Muller  Expansion  for  Symmetric  Functions” 

Jon  T.  Butler,  Gerhard  W.  Dueck,  Vlad  P.  Shmerko,  and 
Svetlana  Yanuskevich 


Abstract — The  above  paper^  finds  an  optimal  fixed-polarity  Reed-Muller 
expansion  of  an  n -variable  totally  symmetric  function  using  an 
OFDD-based  algorithm  that  requires  O(n^)  time  and  0(t»®)  storage 
space.  However,  an  algorithm  based  on  Suprun’s  transeunt  triangles 
[1],  [3],  [4]  requires  only  O(n^)  time  and  O(n^)  storage  space.  An 
implementation  of  this  algorithm  yields  computation  times  lower  by  several 
orders  of  magnitude. 

Index  Terms — FPRM  (fixed  polarity  Reed-Muller  expressions),  two-level 
AND/EXOR  forms,  symmetric  functions,  logic  synthesis,  minimization. 


I.  INTRODUCTION 

A  recent  program,  Sympathy,'  for  finding  optimal  polarity 
Reed-Muller  (FPRM)  expansions  of  symmetric  functions  is  based  on 
an  algorithm  whose  data  structure  is  an  OFDD  of  the  given  function. 
It  requires  0(n‘ )  operations  and  0(n^)  storage  space,  where  n  is 
the  number  of  variables.  However,  if  one  uses  a  more  efficient  data 
structure,  specifically  the  transeunt  triangle  of  Suprun  [1],  [3],  [4], 
the  same  computation  can  be  done  with  0(n^  )  operations  and  0(n^ ) 
storage  space.  The  improvement  is  achieved  because  coefficients 
needed  in  various  expansions  are  computed  and  stored  only  once, 
whereas  Sympathy  builds  a  new  OFDD  for  each  polarity.  On  bench¬ 
mark  functions,  the  speed  improvement  is  by  orders  of  magnitude 

II.  Notation 

A  FPRM  expansion  for  a  general  function  f{xi,  X2,  . . . ,  Xn)  is 

f{xi,  X2,  .  .  .  ,  Xn)  =  CoffiClxt  ®Ci2:2  ffi-  ■  -ffiCnX*  Qc„  +  ixlx2 

S)  ■  ■  ■Q)C2n-lxlx2  ■  ■  ■  Xn  (1) 

where  x*  is  either  Xi  or  Xi  everywhere.  The  term  fixed-polarity  refers 
to  the  fact  that  each  variable  occurs  in  the  expression  in  only  one  way, 
Xi  or  Xi.  For  example,  f{x±,  X2,  xs)  =  x±X2X3  +  x±X2X3  has  the 
following  four  FPRM  expansions. 

No  variables  complemented:  Iffifxi  ®X2  ®X3]®[xiX2  ffixixs  ffi 
2:22:3] 

One  variable  complemented:  xi  ffi  2:12:2  ffi  2112:3  ffi  2:22:3 
Two  variables  complemented:  2:3  ®xiX2  ffi  Ti2:3  ffi  X2X3 
All  variables  complemented:  1  ffi  [2I1  ffi2l2  ©2^3]  ffi  [2;i2;2  ©2:12:3  ffi 
X2X3] 

Note  the  total  number  of  product  terms  required  to  realize  this  func¬ 
tion.  In  the  first  and  fourth  FPRM  expansions,  seven  terms  are  required, 
while  in  the  second  and  third,  only  four  are  required.  The  FPRM  sim¬ 
plification  problem  is  to  determine  which  of  n  + 1  polarities  (number  of 
complemented  variables)  yields  the  FPRM  expansion  with  the  fewest 
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Fig.  1.  The  transeunt  triangle  for /(xi,  X2,  X3)  =  X1X2X3  -|-  X1X2X3. 
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Fig.  2.  Reed-Muller  expansion  matrices  embedded  in  the  transeunt  triangle  of 

f(xi,  X2,  X3)  =  X1X2X3  -F  X1X2X3. 


terms.  In  this  example,  the  two  middle  polarities  are  both  optimum, 
yielding  an  expansion  of  four  terms  each. 

A  function  f(xi,  2:2,  ■■■x„)  is  (totally)  symmetric  if  and  only 
if  it  is  unchanged  by  any  permutation  of  variables.  For  example, 
f(xi,  X2,  X3)  =  X1X2X3  +  X1X2X2  is  symmetric.  Certain  coeffi¬ 
cients  in  the  FPRM  expansion  of  a  symmetric  function  are  identical. 
Let  the  Reed-Muller  expansion  matrix  of  a  symmetric  function  be  an 
(n  +  1)  X  (n  +  1)  matrix  of  binary  coefficients 

doo  dor  •  •  •  don 
dio  dll  ■  ■  ■  din 


BMi  = 


(2) 


\_(lnQ  ^nl  ‘  ‘  ‘  ^nn  J 

where  djk  is  the  coefficient  of  a  product  term  of  2:i’s  in  an  FPRM 
expansion  (1)  in  which  j  variables  are  complemented  and  k  are  not. 
For  the  four  FPRM  expansions  of /(xi,  X2,  2:3)  =  ©2  T3 +2:12:22:3, 

we  have 
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The  bold  values  represent  coefficients  in  the  corresponding  FPRM  ex¬ 
pansion.  The  Os  not  in  bold  are  0  in  all  Reed-Muller  expansion  matrices 
for  the  same  polarity. 

A  symmetric  function  is  completely  specified  by  a  carry  vector  of 
logic  values  A  =  [ao,  ai,  . . . ,  a„],  such  that  f{xi,  2:2,  . . . ,  Xn)  is 
Oi  for  all  assignments  of  values  to  2:1,  2:2,  Xn  that  have  i  Is, 
where  0  <  i  <  n.  For  example,  the  carry  vector  of  f{xi,  2:2,  2:3)  = 
2:12:22:3  +  2:12:22:3  is  [1,  0,  0,  1], 


III.  Transeunt  Triangle  Representation  of  Reed-Muller 
Expansions 

Consider  a  triangle  of  Os  and  Is,  where  the  base  is  a  symmetric  func¬ 
tion’s  n  +  1-bit  carry  vector.  Immediately  below  this  is  a  vector  of  n  Is 
and  Os  formed  by  the  exclusive  OR  of  adjacent  bits  in  the  carry  vector. 
Immediately  below  this  is  a  vector  of  n  —  1  Is  and  Os  formed  by  the 
exclusive  OR  of  adjacent  bits  in  the  previous  vector,  etc..  At  the  bottom 
is  a  single  1  or  0.  Doing  this  for  f(xi,  2:2,  2:3)  =  X1X2X3  +  2:12:22:3 
yields  the  triangle  in  Fig.  1. 

The  resulting  triangle  is  the  transeunt  triangle,  originated  by  Suprun 
[3],  [4].  Notice  that  the  bits  along  the  triangle’s  left  side  are  coefficients 
in  R.Mo ,  while  bits  along  the  right  side  are  coefficients  in  RM3 .  Ad¬ 
ditionally,  embedded  rectangles  represent  the  coefficients  in  RM 1  and 
R,M2  .  This  can  be  seen  in  Fig.  2. 
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TABLE  I 

Execution  Times  (in  secs.)  for  Symmetric  Benchmark  Functions* 


Function 

In 

Out# 

Car.  Vec. 

Opt.  Pol. 

Products 

FDD 

Sympathy 

Symphony 

col4 

14 

0 

oio*-* 

14 

14 

368.1 

0.3 

0.000118 

ml 

6 

6 

10® 

6 

1 

1.0 

0.1 

0.000016 

m4 

6t 

8 

10® 

6 

1 

4.2 

0.1 

0.000016 

misex2 

5t 

5 

0*1 

0 

1 

t 

0.4 

0.000011 

misj 

lot 

10 

liOO 

0 

2 

X 

0.6 

0.000050 

xor5 

5 

0 

(01)® 

0,2,4 

5 

0.5 

0.1 

0.000011 

rd53 

5 

0 

on® 

0 

5 

- 

- 

0.000012 

rd53 

5 

1 

0®1®0® 

0,5 

10 

- 

- 

0.000011 

rd53 

5 

2 

(01)3 

0,2,4 

5 

- 

- 

0.000011 

rd73 

7 

0 

04i4 

0 

35 

- 

- 

0.000022 

rd73 

7 

1 

(0®!®)® 

0 

21 

. 

- 

0.000023 

rd73 

7 

2 

(Oil-* 

0, 2,4,6 

7 

. 

- 

0.000022 

rd84§ 

8 

0 

(0®1®)®0 

0 

28 

4.7 

0.1 

0.000029 

rd84§ 

8 

1 

(01)^0 

0,2,4, 6,8 

8 

- 

- 

0.000030 

rd84§ 

8 

2 

0«1 

0 

1 

- 

- 

0.000030 

rd84§ 

8 

3 

O^l'^O* 

0 

70 

5.3 

0.1 

0.000029 

sym4 

4 

0 

01®0® 

4 

6 

- 

- 

0.000008 

sym6 

6 

0 

0®l3o® 

0,6 

36 

- 

- 

0.000016 

sym9 

9 

0 

03j^4g3 

4,5 

173 

11.7 

0.1 

0.000039 

symlO 

10 

0 

0*1*0® 

0 

266 

27.6 

0.2 

0.000050 

3yml2 

12 

0 

0*1*0* 

0,12 

1288 

• 

- 

0.000079 

symlS 

15 

0 

0*1®0* 

7,8 

15139 

0.000139 

dbniijn^ 

4 

0 

0®1®0 

0 

6 

- 

0.000008 

dbruijn.3 

9 

0 

03l0l30® 

7 

256 

- 

- 

0.000039 

dbruijn-4 

18 

0 

0* 10® 1® 0101*0® 

17 

106,284 

- 

• 

0.000297 

dbruijn^ 

35 

0 

0®10®1®0®1010®1®0101®01*0* 

6 

15,215,790,080 

- 

- 

0.002739 

IV.  THE  Algorithm  and  its  Time  and  Space  Complexity 
A.  The  Algorithm 

Note  that  a  single  element  of  the  transeunt  triangle  represents  one 
or  more  coefficients  in  the  various  Reed-Muller  expansion  matrices. 
The  efficiency  of  the  transeunt  triangle  is  due  to  the  fact  that  it  is  not 
necessary  to  recompute  this  coefficient  for  each  polarity. 


Algorithm  1 

1)  Generate 

2)  For  each 
(djk)  f  and 


[4] 

the  transeunt  triangle. 

RMi,  extract  the  coefficients 
compute  the  number  of  product 


terms . 

3)  Choose  an  RM^  with  the  fewest  product 
terms . 


B.  Time  and  Space  Complexity 

The  following  lemma  gives  both  the  time  and  space  complexity  of 
the  above  algorithm.  The  time  complexity  is  due  to  [4]. 

Lemma  4.1:  Algorithm  1  is  an  0(n^)-time  algorithm  that  re¬ 
quires  O(n^)  storage  space  for  computing  the  optimal  fixed-polarity 
Reed-Muller  expansion  of  a  symmetric  function  on  n  variables. 

Proof:  In  applying  the  algorithm,  0(n^ )  storage  locations  are  re¬ 
quired  for  the  coefficients  in  the  triangles.  0(n)  locations  are  required 
to  store  the  number  of  product  terms,  one  for  each  of  the  n  -|-  1  polari¬ 
ties,  for  a  total  of  0(n^ )  locations.  ■ 

The  OFDD  approach  has  time  complexity  0(n‘ )  and  space  com¬ 
plexity  0(n®).  Thus,  Algorithm  1  represents  a  significant  improvement. 

V.  Experimental  Results 

Suprun  [3],  [4]  did  not  apply  his  algorithm  to  benchmark  functions. 
Our  implementation  is  called  Symphony,  (symmetric  phunction 
optimizing  system),  which  is  written  in  C-|— I-  and  compiled  under 


Microsoft’s  Visual  Studio  Version  6.0  for  Windows98.  It  was  run  on  a 
400  MHz.  Pentium  system. 

A.  Comparison  of  Symphony  on  benchmark  functions 

Table  I  shows,  for  certain  symmetric  benchmark  functions,  the  exe¬ 
cution  time  of  Symphony  compared  to  Sympathy  and  to  FDD,  another 
OFDD-based  minimizer  that  does  not  consider  symmetry  [2].  Table  I 
also  shows  the  number  of  inputs  (In),  the  Output  Number  (Out),  the 
Carrier  Vector  expressed  as  a  regular  expression  (Car.  Vec.),  the  po- 
larity(ies)  that  produced  the  optimal  realization  (Opt.  Pol.),  and  the 
number  of  product  terms  in  the  optimal  solution  (Products).  The  three 
execution  times  (FDD,  Sympathy,  and  Symphony)  are  shown  in  sec¬ 
onds. 

As  can  be  seen.  Symphony  is  very  fast,  requiring  no  more  than 
0.0002  secs,  on  any  of  the  functions  considered  by  Dreschler  and 
Becker.  Indeed,  these  execution  times  are  less  than  the  time  interval 
between  real  time  clock  interrupts.  As  a  result,  timing  functions  in 
C-|— I-  return  zero  elapsed  time  for  program  execution.  To  achieve  the 
necessary  resolution,  each  function  was  minimized  2000000  times 
and  the  total  time  was  divided  by  2  000  000. 

Each  dbruijn_k  entry  in  Table  I  is  a  d’Brujin  sequence  indexed  by 
k.  That  is,  each  sequence  contains  exactly  one  copy  of  each  of  the  2^ 
binary  fe-tuples.  Overall,  it  contains  a  total  of  2*^  -I-  fe  —  1  bits.  This 
sequence  is  such  that  decision  diagram  representations  for  such  func¬ 
tions  will  have  many  nodes,  as  there  are  few  repeated  subsequences. 
As  a  result,  algorithms  based  on  decision  diagrams  will  require  more 
computation  time  than  for  other  symmetric  functons. 

Table  II  shows,  for  certain  symmetric  functions  that  are  also 
threshold  functions,  the  relative  execution  times  of  EDD,  Sympathy, 
and  Symphony.  Again,  Symphony  is  fast. 

VI.  Conclusion 

Rather  than  computing  the  entire  EPRM  expansion  for  each  polarity. 
Symphony  computes  and  stores  expansion  coefficients  only  once. 


1388 


IEEE  TRANSACTIONS  ON  COMPUTER-AIDED  DESIGN  OE  INTEGRATED  CIRCUITS  AND  SYSTEMS,  VOL.  19,  NO.  1 1,  NOVEMBER  2000 


TABLE  II 

Execution  Times  (in  secs.)  for  Symmetric  Threshold  Functions 


Function 

la 

Out# 

Car.  Vec. 

Opt.  Pol. 

Products 

FDD 

Sympathy 

Symphony 

t/iredo,in 

10 

0 

0 

11 

21.0 

0.1 

0.000090 

threa  10, 11 

11 

0 

0i0i2 

0 

11 

42.1 

0.2 

0.000114 

t/ire8ii,i2 

12 

0 

0lli2 

0 

13 

87.3 

0.2 

0.000143 

t/ireai2,i3 

13 

0 

Oi*!!* 

0 

13 

177.5 

0.2 

0.000171 

t/ireai3,i4 

14 

0 

0l3i2 

0 

15 

364.2 

0.3 

0.000209 

thre8i4,i5 

15 

0 

qU^2 

0 

15 

756.5 

0.3 

0.000251 

t/irC5l6,l6 

16 

0 

0l5i2 

0 

17 

1555.8 

0.3 

0.000295 

t/ireai6,i7 

17 

0 

0l6i2 

0 

17 

3140.3 

0.4 

0.000346 

tftreair.is 

18 

0 

0l7i2 

0 

19 

6483.9 

0.4 

0.000406 

tAre3i8,i9 

19 

0 

0‘8l2 

0 

19 

13033.1 

0.4 

0.000468 

t/ireai9,2o 

20 

0 

0lfll2 

0 

21 

25870.1 

0.4 

0.000540 

thre320,2i 

21 

0 

O20i2 

0 

21 

52549.4 

0.5 

0.000612 

thres9A,^^i 

22 

0 

02112 

0 

23 

t 

0.5 

0.000691 

thres22,23 

23 

0 

022i2 

0 

23 

t 

0.5 

0.000781 

24 

0 

023i2 

0 

25 

t 

0.6 

0.000883 

thres24,2i 

25 

0 

024i2 

0 

25 

t 

0.6 

0.000996 

thre326,26 

26 

0 

0«1* 

0 

27 

t 

0.7 

0.001109 

27 

0 

026i2 

0 

27 

t 

0.7 

0.001230 

threa27,2S 

28 

0 

027i2 

0 

29 

t 

0.7 

0.001343 

thres28,29 

29 

0 

0“ia 

0 

29 

t 

0.8 

0.001492 

thres29,3o 

30 

0 

029i2 

0 

31 

t 

0.8 

0.001658 

using  the  transeunt  triangle,  and  extracts  them,  as  needed,  to  form 
the  various  expansions.  In  this  way,  it  achieves  a  major  savings  in 
computation  time  and  storage  over  Sympathy,  which  computes  a 
decision  diagram  for  each  polarity. 

An  abbreviated  version  of  Symphony  can  be  accessed  at  http://www. 
oc.nps.navy.mil/~butler/transeunt.html  (word  length  restrictions  on  the 
server  preclude  carrier  vectors  with  more  than  3 1  bits).  Users  can  input 
a  carrier  vector  and  see  the  transeunt  triangle  along  with  the  number  of 
product  terms  for  each  polarity. 
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